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ABSTRACT 

The NINA apparatus, on board the Russian satellite Resurs-01 n.4, 
has been in polar orbit since 1998 July 10, at an altitude of 840 km. 
Its main scientific task is to study the galactic, solar and anomalous 
components of cosmic rays in the energy interval 10-200 MeV n _1 . 

In this paper we present a description of the instrument and its basic 
operating modes. Measurements of Galactic Cosmic Ray spectra will also 
be shown. 

Subject headings: Cosmic rays, Isotope Composition, Energy Spectrum, 
Satellite, Silicon Detector 



1. Introduction 

With the launch of the telescope NINA on 1998 July 10 a wide program of 
satellite cosmic ray observations began. The program aims to study the cosmic ray 
radiation in a broad energy spectrum (from 10 to 10 5 MeV n _1 ) using a serious of 
dedicated satellite missions. 

NINA (a New Instrument for Nuclear Analysis) has been developed by a 
joint program of the Italian National Institute of Nuclear Physics (INFN) and 
the Moscow State Engineering and Physics Institute (MEPhI). INFN consists of 
several Italian Institutes and Universities (WiZard group) who have carried out, 
together with European and American partners, balloon-borne experiments for the 
detection of cosmic antiparticles since 1989 (Golden et al. 1994; Golden et al. 1996; 
Hof et al. 1996; Barbiellini et al. 1996; Boezio et al. 1997; Boezio et al. 1999; 
Basini et al. 1999; Boezio et al. 2000). 
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The link with the Russian counterpart was established in 1994, when the two 
sides started a collaboration and conceived the Russian- Italian Missions (RIM), of 
which NINA is the first step. 

NINA's goal is to detect cosmic ray nuclei of galactic, solar and anomalous origin, 
at 1 AU, from hydrogen to iron, between 10 and 200 MeV n" 1 . The experiment is 
carried out on board the satellite Resurs-01 n.4, developed by the Russian space 
company VNIIEM. The spacecraft was launched into a polar sun-synchronous orbit 
of altitude 840 km (Casolino et al. 1999a; Sparvoli et al. 2000). 

NINA has been joined in space by a twin detector (NINA-2), placed again in a 
polar orbit but at a lower altitude (450 km). NINA-2 is housed on board the Italian 
satellite MITA (Casolino et al. 1999b), launched on 2000 July 15 from the Plesetsk 
launch facility in Russia by means of a Cosmos launcher. This second mission is 
intended to last for three years. 

The RIM missions will then continue with the deployment of the PAMELA 
magnet spectrometer, which will be installed onboard the satellite Resurs 01 n.5 and 
put in orbit at the beginning of the year 2003. The main objective of PAMELA 
is to perform high precision measurements of antiparticle spectra (positrons and 
antiprotons) in the energy range from 0.1 GeV up to 200 GeV. In addition, it will 
measure electrons, protons and the nuclear components of cosmic rays, and will 
search for cosmic antinuclei (The Pamela Coll. 1999; Adriani et al. 1997). 

This article reports on the NINA mission, its scientific tasks, the organization 
of the detector with its ancillary instruments, the interface with the satellite, the 
launch phase, and finally the performance of the telescope in flight. In addition, it 
presents the reconstruction of the energy spectrum of the galactic component of 4 He, 
12 C, and 16 in solar quiet conditions. 
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2. Scientific overview 

NINA has been built in order to investigate the nuclear and isotopic composition 
of low energy cosmic particles. Its low-altitude polar orbit (about 1.1 Earth-radii) is 
particularly suitable for performing observations of particles of different origin while 
traversing regions of different geomagnetic cut-off. The Earth's magnetic field is 
utilized as a spectrometer. According to the coordinates along the orbit where the 
particles are detected, it is possible to make inferences about their origin which can 
be galactic, solar, or anomalous. 

Galactic Cosmic Rays 

Galactic Cosmic Rays (GCR) are a directly accessible sample of matter coming 
from outside the Solar System. The GCR energy spectrum can be well represented 
by a power-law energy distribution for energies above 1 GeV n -1 , but at lower energy 
shows a strong attenuation due to the interaction between the Solar Wind and the 
cosmic particles (Wiedenbeck & Greiner 1980). This is one of the reasons why GCR 
investigations below 200 MeV n _1 have been relatively scarce in the past. 

NINA started its mission in a period of medium solar activity; the next solar 
maximum is foreseen for the year 2000. Due to its technical characteristics and 
its good energy, mass and angular resolution, the telescope is particularly suited 
for exploring the low energy component of the cosmic radiation. The detector can 
record GCRs of very low energy (from 10 up to 200 MeV n -1 ) in the polar sectors of 
the orbit, where geomagnetic effects are virtually negligible. 

Solar Energetic Particles 

The most complete measurements of elemental abundances in the solar corona 
come from measurements of high energy particles accelerated in the large Solar 
Energetic Particle (SEP) events. 
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Initially it was thought that the energetic particles in large SEP events were 
accelerated in solar flares. In recent years, however, it has become clear that in 
the large gradual events particles are accelerated at shock waves that are driven 
out from the Sun by Coronal Mass Ejections (CMEs) (Reames 1990; Reames 1993; 
Reames 1995; Gosling 1993; Reames 1998). These shocks accelerate the ions of the 
chemical elements in a fairly equivalent manner. In contrast, particles accelerated in 
impulsive solar flares show specific elemental enhancements produced by resonant 
wave-particle interactions during stochastic acceleration of the ions from the flare 
plasma. 

New spacecraft observations can extend Solar Energetic Particle measurements 
to heavier elements, to rarer elements and to isotopes. NINA can perform SEP 
observations in the polar sectors of the orbit. Its good mass discrimination can help 
in the determination of their composition and therefore in the comprehension of the 
sources and acceleration mechanisms involved. 

Anomalous Cosmic Rays 

Anomalous Cosmic Rays (ACRs) are a low-energy component of interplanetary 
particles that include the elements H, He, C, N, O, and Ar (Klecker 1995; 
Simpson 1995). They are now known to originate from interstellar neutral particles 
that have been swept into the heliosphere, ionized by solar UV or charge exchange 
with the solar wind, convected into the outer heliosphere, and then accelerated 
to energies of ~ 10 MeV n^ 1 or more (Fisk, Kozlovsky & Ramaty 1974). It is 
commonly assumed that the bulk of ACR acceleration takes place at the solar wind 
termination shock (Pesses, Jokipii & Eichler 1981). 

The observation of the anomalous N, O, and Ne ionic charge composition with 
the SAMPEX satellite (Baker et al. 1993) confirmed the theoretical predictions that 
ACRs are only partially charged; more precisely, singly charged ions dominate only 
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at energies below ~ 20 MeV n _1 (Klecker et al. 1998), while at higher energies 
multiply charged ions become more abundant. 

Being only partially ionized, ACRs have a much greater magnetic rigidity (at a 
given energy per nucleon) than either GCRs, which are essentially fully stripped, 
or SEPs, which have charge states characteristic of coronal temperatures. As a 
result, ACRs ca n be observed to much lower invariant latitude with a polar orbiting 
spacecraft like NINA. 

3. The NINA instrument 

NINA consists of the following 4 subsystems: a) the detector (box Dl), 
composed of 32 silicon layers and the electronics for signal processing, b) the 
on-board computer (box D2), a dual microprocessor dedicated to data processing 
and to the selection of the trigger and the acquisition mode configuration, c) the 
interface computer (box E), which rearranges the data coming from box D2 and 
delivers them to the satellite telemetry system, and d) the power supply (box P), 
which distributes the power supply to the different subsystems. 

The weight and electric power of the complete telescope are respectively 40 kg 
and 40 W, in accordance with the constraints imposed by the satellite. To safeguard 
from possible malfunctions and breaks, all electronic systems are global redundant. 

The detector of NINA (box Dl), manufactured by the Italian company Laben, 
is composed by 16 silicon planes. Every plane consists of a pair of n-type silicon 
detectors, 60x60 mm 2 , each read out by 16 strips mounted back to back with 
orthogonal orientation, in order to measure the X and Y coordinates of the particle. 
The strip pitch is 3.6 mm. 
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The thickness of the first pair of detectors, composing the first plane of NINA, 
is (150 ± 15) jum; all the others, instead, are (380 ± 15) fim thick, for a total 
thickness of 11.7 mm of silicon in the whole detector. The indetermination in the 
total silicon thickness comes from the process of manufacturing, and is greater (in 
percentage) for the first two thinner layers. In order to reduce to the minimum the 
thickness of dead area interposed between the silicon detectors, a special ceramic 
(A1 2 3 ) frame, passing only under the lateral strips 1 and 16 and 625 fim thick, 
has been utilized. The role of the ceramic is to sustain the single silicon structures 
and connect them mechanically and electronically, by means of 64 pins, to the 
corresponding mother-boards. A photo of the box Dl is shown in Figure 1. 

Each plane of the detector with its electronics is mounted on an aluminum 
mother-board. The 16 planes are vertically stacked (a sketch of the box Dl is visible 
in Figure 2). The interplanar distance is 1.4 cm except for the first two planes 
which are separated by 8.5 cm in order to improve the incident particles trajectory 
determination. They define the angular aperture of the telescope, which is about 32 
degrees. The 16 planes are modular, so that mechanically and electronically they 
are interchangeable. Below the 16 silicon planes other 4 modules, dedicated to the 
trigger electronics, silicon power supply, analog- digital conversion, and FIFO, are 
placed. 

The 20 plane structure is housed in a cylindrical aluminum vessel of 284 mm 
diameter and 480 mm height, filled up with nitrogen at 1.2 atm. The vessel is 2 mm 
thick, except for a window above the first silicon plane where it is reduced to 300 
jum (Figure 2). The top part of the vessel is rounded, while the bottom part houses 
the connectors for the interfaces with the other parts of the detector. 

The lateral strips (n.l and n.16) of every silicon layer are used for the 
Anticoincidence System (AC); they are read together by the same electronic channel, 
except for those of plane 1 where they are physically disconnected. A total number 
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of 448 electronic channels, out of the 512 available in the box Dl, are used for the 
particle track and energy information; the Anticoincidence System data occupy 
an additional 30 channels, while the remaining 34 are used for housekeeping data 
(16 plane currents, 2 temperatures, 4 voltages, 1 threshold level, 11 ratemeters at 
different depths of the telescope), which monitor the status of the whole instrument. 

The signals produced by the incoming particles in the silicon strips are first 
amplified and shaped. Every plane of the telescope has two 16 channels preamplifiers. 
Data are then converted to the digital format by means of a 12 bit ADC, with a full 
scale of 2800 mip (1 mip being equivalent to 30400 electrons or about 105 keV of 
released energy). The resolution per channel is thus about 0.68 mip ch _1 , equivalent 
to 0.07 MeV ch" 1 . There are two independent lines going to two different ADC's, 
for redundancy reasons. Only one is operating at a given moment. Before the ADC 
there are two gain amplifiers that can be selected depending on the acquisition 
mode: the first provides an amplification of a factor 32 (used only for noise tests) 
and the second of a factor 1 (active for normal acquisitions). 

After conversion by ADC, data (1024 bytes/event) are sent through a FIFO to 
an 8 channel bus interface with the on-board computer (box D2), built again by the 
company Laben. Here all tasks of event processing are performed, before sending 
the data to the interface computer (box E) for mass memory storage. 

The core of the box D2 are two 8086 microprocessors working with a clock 
speed of 4 MHz. In normal conditions both of them are operating in Master-Slave 
mode: the Master microprocessor receives the event from box Dl and performs 
pedestal suppression and data reduction tasks, while the Slave is used to format the 
data, according to the acquisition mode, and send them to the subsystem box E. 
It also selects the trigger logic, implements the Second Level trigger and interfaces 
most of the telecommands with the silicon detector. 

The interface computer (box E) represents the last step of the NINA data 



processing before the records are sent to the satellite for transmission, via telemetry, 
to ground. Two exemplars of box E, for redundancy, have been built, both realized 
by the Russian company VNIIEM. Finally the power supply subsystem (box P), 
made also by VNIIEM, has the function of electrically connecting the satellite with 
its various subsystems. The primary tension comes from the solar panels, and it is 
nominally 27 V (between 24 and 34 V). Starting from this, the box P provides three 
different tensions, two for the analog part (+6V, -6V) and one for the digital (+5V), 
totally independent. 

Further details about the instrument and its performance during a test-beam 
session can be found in (Bakaldin et al. 1997; Bidoli et al. 1999). 

3.1. Operating modes 

NINA can work in different operating conditions, switched automatically or via 
telecommand, which affect the trigger system. In particular: 

1. Two thresholds for the energy deposited in the single silicon layers have been 
implemented: a Low Threshold (L.T.), corresponding to 2.5 mip, and a High 
Threshold (H.T.), corresponding to 25 mip. 

The level of the threshold is fixed by telecommand. As an alternative, a system 
of automatic switching of the threshold, activated by telecommand, switches 
automatically from Low to High Threshold whenever the external rate raises 
above 10 Hz, to prevent the memory being saturated. 

2. The strips 1 and 16 of every silicon layer, except the plane first, are used in the 
Lateral Anticoincidence System. 

The hardware Lateral Anticoincidence can be turned off by telecommand, for 
instance in case of a malfunction of one of the lateral strips. In this condition, 



a software veto system (part of the on-line Second Level Trigger) selects only 
tracks not hitting the lateral strips. This procedure ensures that the Lateral 
Anticoincidence rejection is always effective. 

3. The planes 15 and 16 can be used as Bottom Anticoincidence. The default 
operating mode adopts plane 16 but, in case of need, plane 15 can be selected 
by telecommand. 

The Bottom Anticoincidence can be totally removed by telecommand, allowing 
the detection of particles crossing the whole apparatus. 

The main trigger of the acquisition system is the following: 

TRG Ml = D lx x D ly x ((D 2x + D 2y ) + (D 3x + D 3y )) , 

where denotes a signal above-threshold coming from plane i, along view j 
(j=x,y). The logic OR of planes 2 and 3 provides redundancy in case of a failure of 
plane 2. 

In the default operating mode, this trigger is used together with the Lateral and 
Bottom Anticoincidence ON, in order to ensure the complete containment of the 
particle inside the detector. This is the condition which allows the best energy and 
nuclear discrimination to be obtained by NINA. Moreover, TRG Ml can be used 
with Low or High Threshold, defining two different intervals of nuclei which can 
be detected. In particular, TRG Ml in High Threshold mode removes most of the 
protons from the trigger. This is the most frequent configuration adopted in orbit. 

It is possible to switch, via telecommand, to a second trigger: 

TRG Ml = [D 2x + D 2y ) x (D 3x + D 3y ) x (D 4x + D Ay ) x {D 5x + D 5y ) , 

which is used again in its basic operating mode with the Lateral and Bottom AC 
ON. This trigger, used for particular data taking demands or in case of failure of 



the first plane, increases the acceptance angle, at the expenses of a slight worsening 
of the angular resolution. The combination of TRG M2 and High Threshold again 
excludes most of the protons from the trigger. 

The acceptance window of particles with TRG Ml in full containment regime is 
shown in Table 1, for Low Threshold. The spectrum of nuclei extends from hydrogen 
to iron in the energy interval 10-200 MeV n _1 . 

The flux of particles changes notably along the orbit. A limit to the acquisition 
capability of the instrument has been provided by organizing a system which, in 
high rate conditions, enables the detector to register events in less detail. Every 60 
seconds the processor in box D2 calculates the rate of particles reaching the detector, 
and selects one of the following acquisition modes: 

1. Full- Format mode (counting rate up to 10 Hz). This mode, in which 
the whole event topology is recorded, is the normal working configuration 
outside the Earth's Radiation Belts and in particular out of the South Atlantic 
Anomaly (SAA). It allows the measurement of the energy released by the 
particle in each silicon detector and the storage of this information. 

2. Ei-E tot mode (counting rate > 10 Hz). At high fluxes it is necessary to 
make an optimal use of the mass memory. In this acquisition mode, a Second 
Level trigger, driven by the processor in D2, restricts the event acceptance 
and calculates the energy E\ released in the first plane and the energy E to t 
deposited in the whole detector by the crossing particle. Instead of the whole 
event topology, only E 1 and E tot are stored. 
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4. Orbit operations 

4.1. The satellite Resurs-01 

The class of spacecraft Resurs-01 is designed for meteorological observations, 
investigations of Earth natural resources, and large and small scale ecological 
monitoring of the terrestrial environment. Besides these tasks, the spacecrafts are 
also utilized for the separation and insertion into orbit of small piggy-back space 
vehicles. 

The Resurs-01 is inserted into an almost circular sun-synchronous orbit, having 
the parameters shown in Table 2. The launch of this class of vehicles takes place 
from the Baikonur launch facility in Kazakhstan, by means of a Zenit launcher. The 
total mass of the spacecraft at launch, including the separable micro-satellites, is 
3200 kg. The mass of the payload alone is 1000 kg. 

The service system incorporates a power supply (providing a voltage working 
range from 24 to 34 V), an attitude control and stabilization system (which maintain 
the attitude of the spacecraft in three axis), a command radiolink (designed to control 
the onboard instrumentation by radiocommands), a program-time device (which 
allows the spacecraft to be in autonomous operation for 3 days), a radiotelemetry 
system, and an onboard time and frequency standard (which provides stable high 
frequency and synchronous signals, and a time mark). 

The instrument NINA is housed into Resurs-01 version n. 4 as shown in Figure 
3. The box Dl is mounted on the top side external to the satellite, in such a way to 
point always to the zenith during the flight. The other boxes are located inside the 
body of the satellite. 

The box Dl has two external sensors to measure its temperature. A heater 
keeps it thermal stable when switched off. The temperatures of E and P, and the 
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stability of the satellite power supply, are monitored on board Resurs. All these data 
are merged with the satellite telemetry and sent to the ground. 

4.2. Control in orbit 

The interaction between the ground stations and NINA during operations is driven 
by telecommands, which give the possibility to activate a total of 24 commands. 
Some are dedicated to operations like power switching (ON/OFF), data transferring, 
memory cleaning, and selection of single or dual microprocessor model. The others 
act on the trigger logic or on the acquisition model, as illustrated previously. The 
transmission of specific telecommands can be performed when the satellite passes 
over the ground stations. A response packet with the status of the telecommand 
settings is sent to Earth each time a change in the command buffer has occurred. 

The default telecommand set initializes the acquisition with TRG Ml, automatic 
switching between Low and High Threshold enabled, lateral strips and bottom 
plane in anticoincidence, and automatically switching of the acquisition mode from 
Full-Format to E\-E tot , according to the counting rate. However, it is possible to 
set any combination of trigger logic, threshold level, anticoincidence and acquisition 
mode by telecommand. After two months of operation, we switched the acquisition 
to High Threshold mode, in order to focus our analysis on high Z particles. 

A special onboard device allows pre-programmed combinations of telecommands, 
acting automatically in specific points of the orbit. These combinations permit, 
for instance, electronic calibration procedures to be performed over the equatorial 
regions (where the counting rate is low), or to stop the data acquisition in sectors of 
the orbit with very high counting rates. 

The average volume of data that NINA transfers from the satellite to ground is 
2 MB day^ 1 , corresponding to more than 20000 events. NINA has a 16 MB mass 
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memory available in the onboard memory storage. Since the average mass memory 
occupation in solar quiet periods and outside the South Atlantic Anomaly is around 
1.5 MB day -1 , there is the possibility to accumulate data for a few days or during 
solar events, for a subsequent transmission. 

5. Detector performance in orbit 

The launch of NINA took place on 1998 July 10. The transmission of its 
scientific data started on August 31, after an initial period needed for stabilization 
of the orbit and overall checks of the satellite functionality. The analysis of the 
first sample of data received showed that the instrument performanced well and 
confirmed the functionality of the whole system. 

During one orbit the satellite has a day-night cycle according to its position 
with respect to the Earth and the Sun. Two of the 34 housekeeping data available 
on NINA give internal temperatures sampled at two different heights inside the 
detector. We measured the behavior of the two temperatures in orbit. The excursion 
of their values between light and shadow was less than 1 degree, as required in the 
construction phase. 

Important information about the status of the detector are provided by the 
ratemeters, which are also part of the housekeeping sector. These are indicators of 
the particle flux impinging on different planes of box Dl, which is a function of the 
orbit of the satellite. Low and high flux ratemeters are implemented at different 
heights inside the telescope. In case of intense flux, the high ratemeters section 
provides information while the low ratemeters may saturate. 

Figure 4 shows the behavior of the low ratemeter implemented on plane n.6 
during one typical orbit. The counting rate is given in hertz, with the saturation 
value at about 420 Hz. From the picture one can easily follow the path of the 
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satellite through the different regions of the Earth's magnetosphere. In particular, 
it can be seen how the flux increases at the Poles with respect to the Equator, 
because at high latitudes the terrestrial magnetic field does not effectively prevent 
low energy particles from approaching the Earth. The spikes visible near the Poles 
are due to low energy electrons which fill the Outer Radiation Belt. In the South 
Atlantic Anomaly the magnetic field has a local minimum and thus the low energy 
proton flux reaches very high levels. This is clearly evident by the saturation of the 
ratemeter counter. 

We have examined the stability of some important parameters of the detector 
with time, during the first 6 months of life of the detector. The pedestal values 
remained stable within 1 ADC channel from October 1998 until March 1999. The 
same holds for the voltages, the threshold values, the temperatures. 



6. Galactic Cosmic Ray measurements 

As mentioned before, after the first two months of operation NINA's activity 
was focused on the detection of particles heavier than hydrogen. In section 6.1 we 
discuss the track selection algorithm for Z>1 particles, that we utilized to calculate 
the GCR fluxes of 4 He, 12 C, and 16 0, which are shown in section 6.3. The algorithm 
of isotope identification and the performance of mass discrimination of NINA in 
orbit are presented in section 6.2. 



6.1. The track selection algorithm 

The optimal performance of NINA in terms of charge, mass and energy 
determination is achieved by requesting the full containment of a particle inside the 
detector, using the Lateral and Bottom Anticoincidence System as a veto. In order 
to reject upward moving particles, tracks accompanied by nuclear interactions, and 



events consisting of two and more tracks, an off-line track selection algorithm for the 
data analysis is needed. 

The selection algorithm for nuclei with Z>1, implemented for the analysis of 
NINA flight data, applies six rejection criteria: 

1. Real particles moving downwards and stopping inside the detector have an 
energy deposition that increases along the track. It is natural, therefore, to 
request tracks to deposit in each view a quantity of energy greater than in the 
previous one multiplied by a constant K±, which takes into account the energy 
fluctuations. If 

E(i) < Kt x E(i - 1), 

for any i in the range from the second crossed view to the one with the 
maximum deposit of energy, the event is rejected. 

2. In order to clean the data sample from particles with nuclear interactions, two 
energies for each track are calculated: 

- ^track = sum of the energies released by the particle from the first hit view to 
the view following the one with the maximum deposit of energy; 

- ^residual = total amount of energy left in all the remaining layers. 
The two energies are compared and events with 

^residual ' > -^-2 X Efrack j 

where K 2 is a parameter to be optimized, are rejected. 

3. Double tracks are eliminated estimating two energies for each crossed view % 
along the particle path: 

- E duster (i) = sum of energy released in the strip with the maximum deposit of 
energy and in the two nearest strips; 



- E noise (i) = sum of the energy released in the other strips of the silicon layer. 
If 

Enoiseij) > K 3 X E duster (i) : 

for any of the i crossed views and where K 3 has to be fixed, the event is 
rejected. 

4. Events with the maximum deposit of energy in the first view are rejected. This 
criterion, together with the condition n. 1, selects downwards moving particles. 

5. Events where the maximum energy release in the X view and in the Y view 
are not in the same or between consecutive detector planes are rejected. This 
request helps filtering double tracks. 

6. In order to reduce the number of particles which leave the detector through the 
space between planes, events which release the maximum of the energy deposit 
per silicon layer in the strips 2 or 15, for any of the crossed layers, are rejected. 

To apply this algorithm to the data collected in orbit it was necessary to choose 
the Ki, K 2 , and K 3 coefficients in such a way to efficiently clean the data sample 
from the background, minimizing at the same time the number of good events 
rejected. 

In order to optimize these values we utilized samples of different types of 
particles, obtained from a beam test session at GSI in 1997 (Bidoli et al. 1999), as 
well as flight data. 

If we define as efficiency e the value e = (1 - ^. rejected good events ^ th begt 

J v N. good events ' ' 

optimization of the K 1: K 2 , and K 3 values that we achieved (K 1 = 0.7, K 2 = 0.01 
and K 3 = 0.01) determined an efficiency equal to e = 0.975 ± 0.003 for all particles 
with Z>1. 
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6.2. Isotope identification 



Charge and mass identification procedures may be applied to the events which 
survive the track selection algorithm. 

The mass M and the charge Z of the particles are calculated in parallel by two 
methods, in order to have a more precise particle recognition: 

a) the method of the residual range (Baker et al. 1993; Hasebe et al. 1993; 
Sparvoli et al. 1997; Bidoli et al. 1999). 

In this method, the charge Z is estimated by means of the product E 1 x E tot . 
Here E x is the energy released by the particles in the first silicon detector (two 
layers) of the tower NINA, and E to t is the total energy released in the whole 
instrument. Figure 5 shows the E\ vs E tot curves of particles resulting from 
the fragmentation of 12 C by means of a polyethylene target, obtained during 
a beam test of NINA (Bidoli et al. 1999). The nuclear families lie on different 
hyperbolas E x x E tot = k(Z 2 ). 

Once the charge Z is identified by its E 1 — E tot hyperbole, the mass of the 
particle is evaluated by applying the following formula: 



where AE is the energy lost by the particle in a thickness Ax measured 
starting from the first plane, the parameter a is a constant depending on the 
medium and b has a value between 1.5 and 1.8 in NINA's energy range. A 
precise evaluation of such parameters for each atomic species has been obtained 
both from real and simulated data. 




(1) 



b) the method of the approximation to the Bethe-Bloch theoretical curve. 

With this second method we estimate the mass M and charge Z of the particle 
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by minimizing the following % 2 quantity : 

N 

X 2 = E ( W * ( AE * Cal - AE * ke ° r ) ) 2 ' ( 2 ) 
i=i 

where AEl eal is the energy released by the particle in the i-th view, AEf heor 
is the corresponding expected value, Wi is the weight for every difference 
Wi = A greai , and the sum is extended to the N silicon layers activated by the 

i 

particle, excluding the last one where the particle stops and the fluctuations of 
the energy deposits are generally very big. 

In order to build such a function, it is necessary to follow step by step the 
particle's path, calculating the scattering angles at every layer. This method 
takes into account also the energy losses in dead layers, thus preventing 
systematic shifts on the reconstructed masses. 

For a complete rejection of the background, only particles with the same final 
identification given by the two methods are selected. Finally, a cross-check between 
the real range of the particle in the detector and the expected value according to 
simulation is a consistency test for the event. 

In Figure 6 the reconstructed masses using eq. 1 for helium isotopes detected 
in orbit are compared to the ones obtained from the test-beam data; the sample 
of particles in orbit has been selected during passages over the polar caps, and in 
period of quiet solar activity. The picture shows that the mass resolution of NINA 
in flight is in good agreement with the measurements performed at GSI. 

6.3. Determination of Fluxes 

The analysis presented in this section refers to particles registered by NINA 
in the solar quiet period December 1998-March 1999, detected in High Threshold 



mode and TRG Ml acquisition. In order to select a sample of pure low energy 
(E>10 MeV/n) primary cosmic rays, and avoid the distortions induced by the Earth 
magnetic field, only particles registered at a value of L-shell>6 (L geomagnetic shell) 
were chosen. 

In order to estimate the cosmic ray fluxes it is also necessary to know the 
geometric factor of the instrument as a function of energy, and the exposure time 
in orbit. A correction factor can then be applied to account for energy-loss in the 
aluminum window. 

The geometric factor of the instrument was calculated by means of Monte Carlo 
simulations based on the CERN-GEANT code (Brun et al. 1994). Each simulated 
track underwent the trigger conditions (TRG Ml or TRG M2), with the Lateral and 
Bottom Anticoincidences ON, both for Low and High Threshold mode. Figure 7 
presents the geometric factor G of NINA for 4 He, 12 C, and 16 in High Threshold 
mode over the energy intervals defined by the trigger and the selection conditions 
explained before. 

The incoming energy of the particles was reconstructed by an iterative algorithm 
which is based on the Bethe-Bloch formula. The algorithm works this way: as a first 
step, the total energy E tot is taken as initial energy E in of the particle; with this 
value of initial energy, all energies deposited in every layer along the particle track 
(namely the aluminum window, the silicon detectors and the inter-gap volumes of 
nitrogen) are calculated, and the expected value of the total energy deposited in the 
silicon E^f estimated. 

The value of E^P is then compared with E tot . If their difference is greater than 
0.1 MeV we define a new initial energy as Ej n = E^f + E step , where E step is an 
incremental step energy fixed according to the precision that we want to reach, and 
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perform a new iteration. When finally the condition 

E tot -EZ P < 0.1 MeV 
is fulfilled, the algorithm stops and the initial energy of the particle is identified. 

The differential energy spectra were then determined by the following formula: 

AN(E) 



Flux(E) = 



TeG(E) AE' 



where AN(E) is the number of detected particles with energy between E and 
E + AE, T is the exposure time in orbit for the period under consideration, e is 
efficiency of track selection discussed above (0.975 for nuclei with Z>1), G(E) is 
the average value of the geometrical factor between E and E + AE, and AE is the 
energy bin chosen to plot the flux. 

Figures 8, 9, and 10 present respectively the differential energy spectra for 4 He, 
12 C, and 16 0, measured by NINA in the solar quiet period December 1998-March 
1999. Errors due to energy resolution are negligible since the energy resolution for 
NINA is better than 1 MeV n _1 , which is much less than the width of the energy 
bins in the flux plots. 

In Figure 8 NINA flux of 4 He is plotted together with the results from the 
mission SIS on ACE, about in the same period of observation. Data from SIS on 
ACE (ACE Home Page) belong to a cycle of 27 days from 1999 February 6 to 1999 
March 4, and are the sum of 3 He and 4 He; errors are statistic plus systematic. 

There is a general agreement among the two sets of results. However there are 
differences between the results of NINA and SIS, which can be attributed to the 
different time period (since the flux of helium is known to change significantly over 
the months) and to the fact that the SIS flux include also 3 He. 

In Figures 9 and 10 the NINA differential energy spectra of respectively 12 C and 
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16 are plotted together with results from the missions SIS and CRIS on ACE, both 
referred to the period 1999 February 6 to 1999 March 4 (ACE Home Page). The 
fluxes measured on board NINA are in very good agreement with the ones on ACE. 

7. Conclusions 

The space telescope NINA, launched in 1998, is a silicon detector devoted to 
the study of cosmic rays of galactic, solar and anomalous origin in the energy range 
10-200 MeV n -1 at 1 AU. It is capable of nuclear identification up to iron and 
isotopic discrimination up to nitrogen, allowing important space physics issues, such 
as the composition and energy spectra of cosmic ray particles, to be addressed. 

The first months of data analysis confirmed that the instrument is working 
properly in space; the overall performances of the detector are in good agreement 
with expectations and, in particular, the mass resolution capability reached by NINA 
in space reproduces the one obtained in a beam test session. 

The energy spectra of galactic 4 He, 12 C, and 16 measured by NINA at 1 AU 
have been determined. The analysis of the galactic ratio 3 He/ 4 He together with the 
abundance ratio of the isotopes of hydrogen is in progress, as well as the study of 
the composition and energy spectra of SEP and ACR particles. 

The measurements performed by NINA are important in view of the 
second mission NINA-2 which will complement the observations of NINA 
extending its lifetime to cover a complete solar cycle. The addition of PAMELA 
(The Pamela Coll. 1999; Adriani et al. 1997) will allow the extension of cosmic ray 
observations to energies greater than 200 GeV. 

We acknowledge the Laben company (Italy) for the realization of Dl and D2 
parts of the detector as well as VNIIEM (Russia) for parts E, P and especially for 



the assembly and space operations on the Resurs satellite. 

Finally we acknowledge the Russian Foundation of Base Research, grant 
99-02-16274, who partially supported the Russian Institutions for this work. 
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Table 1: Energy windows for the most abundant particles fully contained in the 
detector NINA, with TRG Ml and acquisition with Low Threshold. 

Particle Z Energy window (MeV n _1 ) 





1 


10 - 48 


4 He 


2 


9 - 50 


7 Li 


3 


11 - 54 


9 Be 


4 


13 - 65 


n B 


5 


15 - 75 


12 C 


6 


17 - 90 


14 N 


7 


18 - 95 


16Q 


8 


20 -103 


19p 


9 


21 - 107 


20 Ne 


10 


23 -117 


28 Si 


14 


28 - 142 


40 Ca 


20 


39 - 175 


56p e 


26 


58 - 195 
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Table 2: Parameters of the Resurs-01 n.4 satellite. 



Characteristic 


Value 


Orbit inclination (deg) 


98.75 


Orbit period (min) 


101.31 


Eccentricity 


1.12xl0~ 3 


3 axis stabilization accuracy (deg) 


1.0 


Average orbit altitude (km) 


840 




Fig. 1. — Photograph of the internal structure of box Dl. 



- 30 - 




i 1 

UJUillLL, ui-LijiLu-LiiJIiL 
P^i fii>i.: = 


Electronic 
Reed Out 
Boards 




"niffll 

HHJ 
ii iiiiin niii 









DATA 



POVER 



[ [Ml | [S| GAS rMl i [ |S| | Hrf -ft 




=1 E=3 



Fig. 2.- 



Sketch of the internal structure of box Dl. 
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Fig. 4. — Low ratemeter counting rate (at plane 6) as a function of time. The two 
black vertical arrows define 1 satellite orbit. 
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Etot (MeV) 



Fig. 5. — Distribution of the energy released in the first plane (Ei) and the total 
energy (E tot ) detected for particles produced in the fragmentation of 12 C at GSI 
(1997). 
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Fig. 6. — Mass distribution, as given by eq. 1, of a sample of helium isotopes collected 
by NINA in orbit (top) and during the beam test session at GSI (bottom). 
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Fig. 8. — Differential energy spectrum for 4 He in the solar quiet period December 
1998-March 1999 measured by NINA, together with data of SIS on board ACE. 
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Fig. 9. — Differential energy spectrum for 12 C in the solar quiet period December 
1998-March 1999 measured by NINA, together with data of SIS and CRIS on board 
ACE. 
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Fig. 10. — Differential energy spectrum for 16 in the solar quiet period December 
1998-March 1999 measured by NINA, together with data of SIS and CRIS on board 
ACE. 



